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We examine the chiral corrections to exotic meson masses calculated in lattice 
QCD. In particular, we ask whether the non-linear chiral behavior at small quark 
masses, which has been found in other hadronic systems, could lead to large cor- 
rections to the predictions of exotic meson masses based on linear extrapolations to 
the chiral limit. We find that our present understanding of exotic meson decay dy- 
namics suggests that open channels may not make a significant contribution to such 
non-linearities whereas the virtual, closed channels may be important. 



1. Introduction. One of the biggest challenges in hadronic physics is to understand 
the role of gluonic degrees of freedom. Even though there is evidence from high energy 
experiments that gluons contribute significantly to hadron structure, for example to 
the momentum and spin sum rules, there is a pressing need for direct observation 
of gluonic excitations at low energies. In particular, it is expected that glue should 
manifest itself in the meson spectrum and reflect on the nature of confinement. 

In recent years a few candidates for glueballs and hybrid mesons have been re- 
ported. For example, there are two major glueball candidates, the /o(1500) and 
/o(1710), scalar-isoscalar mesons observed in pp annihilation, in central production 
as well as in J/ip (/o(1710)) decays However, as a result of mixing with regular 
qq states, none of these is expected to be a purely gluonic state. This is a general 
problem with glueballs - they have regular quantum numbers and therefore cannot 
be unambiguously identified as purely gluonic excitations. On the other hand, the 
existence of mesons with exotic quantum numbers (combinations of spin, parity and 
charge conjugation, J'^'-^, which cannot be attributed to valence quarks alone), would 
be an explicit proof that gluonic degrees of freedom can indeed play an active role at 
low energies. 

Current estimates based on lattice QCD as well as QCD-based models [Q, 

suggest that the lightest exotic hybrid should have J"^*" = 1 ^ and mass slightly 
below 2 GeV. On the experimental side, two exotic candidates with these quantum 
numbers have recently been reported by the E852 BNL collaboration. One, in the rj-K^ 
channel of the reaction 7r~(18GeV)p Xp —>■ rjir^p, occurs at Mx ~ 1400MeV [^, 
while the other, in the and ri'ir channels, has a mass Mx ~ 1600MeV P, 0. 
Although the lighter state could be obscured by final state interaction effects as well 
as leakage from the strong X = a2 decay, the heavier one seems to have a fairly clean 
signal. 

In the large Nc limit it has been shown that exotic meson widths obey the same 
Nc scaling laws as the regular mesons |^. The lack of an overwhelming number of 
exotic resonance candidates may therefore be the result of a small production cross- 
section in the typical reactions employed in exotic searches - e.g., high energy vrA^ 
or KN scattering. Indeed, low lying exotic mesons are expected to have valence 
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quarks in a spin-1 configuration and their production would therefore be suppressed 
in peripheral pseudoscalar meson-nucleon scattering. By the same arguments one 
would expect exotic production to be enhanced in real photon-nucleon scattering P, 



10, 11]. This is particularly encouraging for the experimental studies of exotic meson 



photoproduction which have recently been proposed in connection with the JLab 
12 GeV energy upgrade |T2[ . 



In view of the theoretical importance of the topic and the exciting new experi- 
mental possibilities for exotic meson searches it is of prime importance to constrain 
the theoretical predictions for exotic meson masses - especially from lattice QCD. 
One of the important practical questions in this regard is the effect of light quarks 
on the predicted spectrum. All calculations that have been performed so far involve 
quarks that are much heavier than the physical u and d quarks and, with the excep- 
tion of the initial study by the SESAM collaboration (based on full QCD with two 
degenerate Wilson fermions), are performed in the quenched approximation. This is 
directly related to the present limitations on computer power. It is well known that, 
as a consequence of dynamical chiral symmetry breaking, hadron properties are non- 
analytic functions of the light quark masses. This non-analyticity can lead to rapid, 
non-linear variations of masses as the light quark masses approach zero. Whether 
or not this could make a sizable difference to the lattice QCD predictions for exotic 
state masses is the question we address. 

There has been considerable activity concerning the appropriate chiral extrap- 
olation of lattice results for hadron properties in the past few years, ranging from 



magnetic moments [T^ to charge radii [14] and structure functions [15] as well as 
masses 



T6| . The overall conclusion from this work is that for current quark masses in 



the region above 50-60 MeV, where most lattice results are available, hadron prop- 
erties are smoothly varying functions of the quark masses - much like a constituent 
quark model. However, as one goes below this range, so that the corresponding pion 
Compton wavelength is larger than the source, one finds rapid, non-linear variations 
which are a direct result of dynamical chiral symmetry breaking. These variations 
can change the mass of hadrons extracted by naive linear extrapolation by a hundred 
MeV or more. (Certainly this was the case for the N and A [|T^, while for the p the 
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difference was only a few MeV 0].) A variation of tliis order of magnitude for tlie 
predictions of exotic meson masses would clearly be phenomenologically important 
and our purpose is to check whether this seems likely. 

2. Chiral corrections from pion loops. Our investigation of the possible non- 
linearity of the extrapolation of exotic meson masses to the chiral limit, as a function 
of quark mass, will follow closely the earlier investigations for the and A baryons 
and the p meson. The essential point is that rapid non-linear variations can only arise 
from coupling to those pion-hadron channels which give rise to the leading (LNA) or 
next-to-leading non-analytic (NLNA) behavior of the exotic particle's self-energy. In a 
decay such as E ttH, this means that H must be degenerate or nearly degenerate 
with the exotic meson E. In addition, the relevant coupling constant should be 
reasonably large. For example, in the case of the p meson the relevant channels are 
UTT (LNA) and nn (NLNA) and the extrapolation of the p meson mass was carried 
out using: 

mp = a + bml + a^^ + a^^, (1) 

where the self-energy contributions from the channels ij = un or rnr. 

Our first step in studying the extrapolation of exotic meson masses is therefore 
to look at the channels to which the exotics can couple which involve a pion and 
to check what is known about the corresponding coupling constants. The matrix 
element describing a decay of the J^^ = 1 isovector, exotic mesons with mass 
mx, is given by: 

LS Ml Ms V I 

xgLs{k{mx)){sAXA,SBXB\SMs){SMs,LML\lMx){lAlllBl'B\^,l'x), (2) 

where k(mx) = |k| is the break-up momentum of the AB meson pair, from a decay 
of a state of mass mx, 



,f ^ ,'(m|-(mA-mij)2)(m|-(mA + mB)2) 

k{mx)=\ ^^^^2 ) = Hmx^mA^niB), (3) 

produced with angular momentum L and spin S, with L + S = 1 
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Decay channel 


wave 


PSS 

gl{k{mx))/An F^l MeV] 


IKP 

gl{k{mx))/A7T r^MeV] 


r]TT 

1] 7T 

pn 
/2(1270)7r 
/i(1285)7r 

6i(1235)7r 

r/„(1295)7r 
p(1450)7r 


P 
P 
P 
D 
S 
D 
S 
D 
P 
P 


9.4 X 10-3 0(10-2) 
9.4 X 10-3 0(10-2) 
8.32 0(10) 
5.3 0(10-2) 
1.6 0(10) 
2.3 0(10-2) 
6.1 O(IOO) 
37.9 0(1) 

37.6 0(10) 

30.7 0(10-2) 


7.3 X 10-3 0(10-2) 
7.3 X 10-3 0(10-2) 
5.95 0(10) 

1.9 0(10) 
309. 0(1) 
6.54 O(IOO) 
324. 0(10) 
21.25 0(10) 
15.8 0(10-2) 


Table 1: Isovector Hybrid decay parameters 



In terms of the couphngs, qls, the partial decay widths are given by 

/ 1 \ 2^+1 



327r2 \mx 



(4) 



These couphngs have been calculated in Refs. |]19|, in two models based on the flux 
tube picture of gluonic excitations but assuming different qq production mechanisms. 
In general the models agree on predicting sizable couplings to the so called "S+P" 
channels, where one of the two mesons in the flnal state has quarks in relative S- 
wave and the other in relative P-wave, e.g., nbi. In Table 1 we summarize the 
results of these calculations, listing only the decay channels containing a pion (S-wave 
quarks). Since the overall normalization of the decay matrix elements in these models 
is somewhat arbitrary, we have rescaled the original predictions given in Ref. to 
match the total width ( 170 MeV) and mass {mx = 1.6 GeV) of the exotic pyr state 
found by the E852 experiment [Q. 

Two features deserve particular attention. The couplings to ground state meson 
multiplets are generally smaller than to excited meson multiplets. This may be ar- 
tificial; it is associated with the underlying structure of the model for pair creation 
which, when applied to exotic decays, strongly suppresses final states with mesons 
which have similar orbital wave functions. These models are based on a simple quark 
model picture in which (say) the vr, p and rj have very similar orbital wave functions. 
This is obviously an oversimplification. Secondly, the PSS and IKP models are dras- 
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Decay channel 


wave 


gl{k{mx))/A'n 


Vl[ MeV] 




P 


24.59 


O(IOO) 


6i(1235)7r 


S 


7.12 


O(IOO) 



Table 2: Isovector Hybrid couplings g\{k{mx)) /^t^ 

tically different when it comes to predicting ratios of branching ratios. In general, 
PSS predicts that higher partial waves should be strongly suppressed. 

An alternative approach was presented in Ref. There, it was assumed that 
the decay of the observed exotic is dominated by just two channels : p7r - the one in 
which it has been seen - and the rest, say hiTi (in an S-wave). The couplings obtained 
this way are shown in Table II. 

The shift in the mass of the exotic meson associated with the pion loop, S, can 
be calculated in second order perturbation theory (from S = (P'|l^G'l^|P')/(P'|P)) 
and is given by: 



E 

L 



mi 



dM 



QA-n^ JmA+rriB M \ ITLx j 
V ^ ^^^^ 



V'"-' gliHM)) 



647r3 



mx - M 
k V 



am 



(5) 



with k{M) = \{M^mAiTnB)- The leading non-analytic (LNA) behavior of this self- 
energy is obtained by extracting the piece, independent of the ultra-violet cut-off (or 
form factor, g^lk) / gL{k{mx))) , which exhibits the strongest non-analytic variation as 
a function of the quark mass as one approaches the chiral limit. (This is the term with 
the lowest odd power of or the lowest power of m,r multiplying Inm.,^.) Setting 
niB = rriTr the form of the LNA behavior depends on the relation between mx — ttla 
and and can be easily extracted analytically from Eq. (^ in two limiting cases. 
Consider first the limit mx — ttla << m^, corresponding to an off-shell transition 
between the exotic meson and a nearly degenerate meson plus pion. In this case the 
self energy contribution reduces to 



E 

L 



-V / dkk'^ I 

647r'^ JO \mx 



k"^ + ml ' 



(6) 
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and the LNA behavior is given by 

hLNA = 2^[-l) — • {') 

^ SzTT 47r \mxJ 

and /ctt ~ 0. In the case mx — >> ^-n-, corresponding to a physical decay process 
to two hght mesons (one of them being the pion), one obtains, 



) mAijnx — mA)\Jk? + 



leading to 



m\ {2L + iy\gl{K) /m!^^+' 



^ IQ'K^mAymx — \2,L + 2)\\ An ym^ J 

Even though, in the present case, neither of these is the exact, four-dimensional, pion 
loop contribution, they should give a reliable guide as to the non-linearity that one 
may expect in the extrapolation of lattice results to small quark mass. 

In order to estimate the full self energy contribution, we need to know the off-shell 



dependence of the coupling constant. In the models of Ref. [Tl|, |T^, the momentum 
dependence of the couplings arises from the Fourier transform of quark wave function 
overlaps and is typically of a gaussian form, 

gL{k) = gL{k{mx))e-'"/"'"+''"^^'^^/^''" (10) 

with the scale parameter expected to be in the range (3 ~ 0.2 — 1 GeV. To be 
consistent with the approximation of a heavy source, corresponding to Eqs. (^) and 
(I), we need to set k'^imx) ~ or k'^imx) ~ {jnx — ttia)'^ respectively. Thus the 
momentum dependence of the form factor near k"^ = k^ is not expected to significantly 
renormalize the value at k"^ = k'^{mx) and in Eqs. (^ and (P) we can simply use the 
on-shell couplings from Tables 1 and 2. However, the off-shell transition matrix 
element may have a more significant momentum dependence. This could happen, 
for example, because of nodes present in the radial wave functions of excited mesons 
{e.g. the p(1450) ). Such nontrivial momentum dependence of the form factor could 
possibly introduce a factor of 2 — 3 uncertainty in our predictions. Keeping this 
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Figure 1: Pion loop contribution to the exotic meson self energy from P-wave 
couplings calculated for (a) /3 = 400 MeV (left) and (b) /3 = 700 MeV (right). The 
contributions from open and closed channels are shown with the dashed and dotted 
lines, respectively and the net correction is shown with the solid line. 

possibility in mind, using Eq. (|^) we find that for the maximal strength P-wave 
decays with couplings of the order of 20, 



S ^ -O.lmf 



with both S and in GeV, and all couplings of order 1 are irrelevant. For the 
D-waves with coupling as large as 300 (IKP model, D-wave in the bin channel) we 
obtain, 

S ^ +0.5m^, (12) 

from Eq.( |7p or S ^ +0.4m^lnm,r from Eq. (|^). In the case of the bin or /ivr decay 
channels the second approximation is probably more accurate. 

To be more precise we also check how the full formula for the energy shift com- 
pares with the simpler expression for a heavy source given above. In Fig. 1 the dashed 
line shows S as a function of the pion mass, calculated including the three largest 
P wave open channels from Table 1 (p(1450)7r, ?7„(1295)7r and pvr) using the PSS 
(larger) couplings. The magnitude of S is determined by the scale parameter cho- 
sen as p = 400 MeV in Fig. 1(a) and p = 700 MeV in Fig 1(b). For small p the 
dominant contribution comes from the lightest open channel. This is because there 
is a large mismatch between k{mx) and (k) ~ /5 which leads to enhancement from 
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Figure 2: Pion loop contribution from the open, D-wave channels. The shaded 
region corresponds to (3 in the range from 400 to 700 MeV 



gL{k) (Eq. ([T0|) ). As (3 increases it is the channels with largest on-shell coupling that 
dominate. In Fig. 2 the contribution to S from the two largest D — wave couplings 
(/ivr and feivr in the IKP model) is shown for (3 in the range between 400 MeV and 
700 GeV. In all cases we find that the mass shifts from the open channels do not 
exceed 50 MeV. The magnitude of the self energy is naturally sensitive to the scale 
parameter /?, in particular for higher partial waves. Nevertheless, the weak m^r depen- 
dence for all (3, in particular for the two limiting values corresponding to the upper 
and lower bounds, decreasing for higher partial waves, is consistent with Eqs. (|1T]) 
and dg). 

These are small shifts as compared, for example, to the case of the pion contri- 
bution to the nucleon mass, where the LNA term is of order — 5.6m^. To be more 
specific, the ttNN vertex, written in the notation of Eq. (^, becomes 

X IX X 

g^NNu{p'^ A, /3)75r^? m(p. A, h) = mN9nNN{-h, 1^3 I 0^3) (o l^^l o -^') — ^ia^l i^^) 

Z Z Z Z TTIn 

(13) 

where 



QttNN 



(14) 



W 3 



3 9\ 



(15) 



Eq. (|^ then gives, 

_ 3 X ingl^j^ /^rw \^ _ _ 

327r An \mx J ~ 32tt m% ""^ 327r ' 

Comparing gl^j^/Air = 127r x 14.4 ^ 540 to the typical P-wave coupling from Table 1 

or 2, it becomes clear why the open channels give small corrections to the exotic 

meson mass. 
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However, just as in the case of the nucleon, one expects the largest LNA behavior 
to come from virtual transitions to nearly degenerate states; in this case from tran- 
sitions between an isovector exotic and a pion plus an isoscalar exotic meson. The 
couplings in Tables 1 and 2 only account for real decays and therefore cannot be used 
to estimate such transitions. To the best of our knowledge there are no microscopic 
calculations of such matrix elements, however, these can in principle be derived from 
PCAC In the soft pion limit one has, 

y, p', / = 1, a\q,A^\m'^. A, p, / = 0) = 

p'. A', / = 1, a| V|7r^ q; 1"+, A, p, J = Q)/{2Txf5\^' - p). (16) 

If one assumes the flux tube does not affect the axial charge, then in the static limit 
the valence quark contribution to the Ihs is given by 2mx5a6e*(A') ■ [e(A) x q\. After 
comparing with Eq. (^ this yields 

5'l-+(/=l)^l-+(7=0),7r = 2y — — . [i I ) 

\ 'i In 

The lack of any contribution from the fiux-tube makes this coupling identical to 
the one for ordinary mesons, e.g. the ptuvr coupling (with rrix = rnp ~ wi^). The 



phenomenological value for the gpujn, in our notation, is gpuj-n/mp = ySn/Sgp^T, with 
gpujTT ~ 15 GeV~^, so that the simple quark model overestimates the coupling by 
approximately 50%. The flux-tube may contribute to the axial current if it couples 
to the spin of the quarks. In particular, since it is expected that the ground state 
corresponds to a fiux-tube in a P-wave with respect to the valence QQ pair, the 
overlap of the spin and orbital wave functions may modify the numerical factor in 



Eq. (|T3)- However, it is not expected to alter the ySir/Smx/fTv enhancement arising 
from the soft pion emission. Thus, from Eqs. (^ and ([T7| ) we estimate 

S = -T^2^l ^ -3^'- (18) 

The magnitude of this correction is almost as large as that found in the nu- 
cleon case. The contribution to S from the virtual transitions calculated using the 
expression in Eq. (B|), is shown in Fig. 1 with the dotted line. Its magnitude is 



9 



strongly dependent on the cutoff parameter jS, and can be as large as (9(100MeV) 
for /3 = 700 MeV. It also has significant variation as a function of expected 
from Eq. (p!8|). The total self energy shift arising from the real and virtual corrections 
is shown in Fig. 1 with the solid line. This simple analysis of virtual corrections 
also applies to the p(1600)7r channel which, depending on the relative phase, could 
enhance the overall LNA behavior by a factor of two. Finally if the other light exotics 
J^^ = O"*"", 2^ are not too far from the J^*" = 1 ^, they could significantly enhance 
the (i-wave (oc m^) behavior (c/. Fig. 2). 

3. Conclusions. We have explored the self-energy corrections to the mass of the 
1 ^ exotic meson which are most likely to introduce some non-linearity in the chiral 
extrapolation of lattice estimates of its mass. As in earlier work involving the and 
A baryons and the p meson, the guiding principle was to find those coupled channels 
involving a pion which yield the leading and next-to-leading non-analytic behavior 
of the self-energy. Using several models for the coupling constants to these decay 
channels available in the literature we find that the effects of the open channels are 
rather too small to lead to any significant non-linearity. 

On the other hand, the virtual corrections can potentially be as large as those 
found for the nucleon, A or p meson. Thus current estimates of the masses of exotic 
mesons, based on linear chiral extrapolation of lattice results obtained at relatively 
large input quark masses, may not be sufficiently accurate and nonlinear behavior sim- 
ilar to that found for the nucleon and A seems likely. A careful chiral extrapolation, 
constrained by new calculations at lower quark mass (which are needed urgently), 
may well lead to a physical mass for the 1 ^ of order 100 MeV or more lower than 
that found by naive linear extrapolation. To estimate the strength of the virtual tran- 
sition we used the soft pion limit theorem and in the future more precise calculations 
of the axial current matrix elements involving exotic mesons should be performed to 
test this approximation. Finally, it will be essential to revisit these conclusions as we 
learn more about the open channel exotic decay modes. In particular, if the existing 
models turn out to significantly underestimate these pion couplings, then one would 
also need to revise their contribution to the lattice extrapolation procedure. 
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